Two-Dimensional Nucleation with Edge and Corner Diffusion 
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The effect of edge and corner diffusions on the morphoiogy and on the density of isiands nucieated 
irreversibly on a fiat substrate surface is studied. Without edge and corner diffusion, islands are 
fractal. As an edge diffusion constant D e increases, islands tend to take a cross shape with four nee- 
dles in the < 10 > direction. Additional corner diffusion with a diffusion constant D c yields square 
islands. When D e is small relative to the surface diffusion constant D 3 , the square corner shows the 
Berg instability to produce hopper growth in the < 11 > direction. The corner diffusion influences 
the island number density n. At a deposition flux F with a small D c , mainly monomers are mobile 
and n oc (F/ Ds) 1 ^ 3 . At large D c , dimers and trimers are also mobile and n oc F 3 ^ 7 D7'^ 21 Dc 4 ^ 21 ■ 
The F dependence is in good agreement to the rate equation analysis, but the dependence on D c 
cannot be explained by the theory. 

PACS numbers: 68.55.-a, 68.43. Jk, 05.70.Ln, 81.10.Aj 



I. INTRODUCTION 

Surface morphology during epitaxial growth depends 
on processes taking place on the substratei 1 ^ 3 ^^' 7 '? 
At very low temperatures as in the case of molecular 
beam epitaxy (MBE) , deposited atoms are adsorbed and 
migrate on the substrate surface until they form nuclei 
or are incorporated into islands, but they never evap- 
orate back to the ambient vapor. When an adsorbed 
atom (adatom, for short) sticks and freezes irreversibly 
to an island, the resulting aggregate takes an irregular 
dendritic foim^^^ called diffusion-limited aggregate 
(DLA)iS It has a characteristic feature of self-similarity, 
and is regarded as a typical example of a fractal. In 
actual experiments, one does not always observe fractal 
but also compact islands^ Compact islands are pos- 
sible if they can dissociatoiS or if atoms at the edge 
of them can diffuse along the peripheryi2i2iiiii2iiiii2iiiiii& 
Since at low temperatures it is difficult for an edge atom 
to detach, the edge diffusion mainly governs the island 
morphology** During the edge diffusion, however, the 
outer corner provides an additional barrier to be sur- 
passed. There, an edge atom has to pass an inter- 
mediate state where the number of neighboring bonds 
is less than that on the straight step edge. The edge 
and the corner diffusion is known to influence island 
morphologyi 2 ! 6 ! 7 ! 8 ! 11 : 12 ! 13 : 14 ! 15 : 16 Without edge and cor- 
ner diffusion, islands are similar to DLA. With the edge 
diffusion they take cross shapes, with the corner diffu- 
sion in addition they become square. At the coalescence, 
these peripheral diffusion governs the shape relaxation^ 
Furthermore, the corner diffusion barrier is found to in- 
duce mound formation for multilayer depositions* 1 ^ 

Diffusion processes are also known to influence the 
number density n of islands nucleated on a singular 
surfacei 20 i 21 i 22 i 23 i 24 i 25 i 26 When the deposition rate F is 
large, the adatom density is high and the nucleation oc- 
curs frequently. Thus the island density n is expected to 
be high. When an isolated adatom migrates on the sur- 
face with a large surface diffusion constant D s , the de- 



posited atom moves around a long distance, until it finds 
nucleation partners or preexisting islands and sticks to 
them. Thus, as the diffusion rate D s increases, the rate 
of forming new nuclei decreases, and the island density n 
on a singular surface diminishes. 

More quantitatively, the island density n is defined as 
the total number density of immobile and undissociable 
clusters with sizes larger than the critical cluster size i*. 
In the classical mean-field nucleation theory, the number 
density of adatom clusters of size i is analyzed by means 
of rate equations When clusters smaller than 

i* can dissociate but only the mononer adatom is mobile 
with the diffusion constant D s , n is found to follow the 
scaling law 2 ^ 



(1) 



If the clusters with sizes i smaller than the critical size 
i* cannot dissociate but arc mobile as a whole with their 
respective diffusion constant Di. the scaling form is ap- 
proximately given as 2 - 4 * 2 *^ 7 - 
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By assuming that islands are point sinks to diffusing 
adatoms and the critical nuclear size is i* = 1, the island 
density scaling n oc (F/Dg) 1 ^ 3 is confirmed by the kinetic 
Monte Carlo (KMC) simulation^ Actually, islands in- 
crease their size during growth, and the scaling relation 
between the density and the diffusion constant might de- 
pend on the island morphology. There arc already many 
studies on the island density by means of KMC sim- 
ulation when clustes can dissociate&iiiSiSSi There are 
also KMC simulations with mobile clustersi 27 i 29 i 30 i 31 We 
study here the effect of edge and corner diffusion pro- 
cesses on the island density and its relation to the island 
morphology, systematically 

In § 2, we explain our simulation scheme. In § 3, the 
effect of peripheral diffusion on the island morphology 
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is studied. On varying the ratio between the surface, 
edge and corner diffusion, morphology varies as fractal, 
needle, and square shapes, as previously observed^ We 
point out the different roles of the edge and corner dif- 
fusions on the island symmetry. The island density is 
discussed in § 4. As long as the corner diffusion is small, 
the density scaling is found to be almost independent of 
the island morphology with an exponent 1/3, in consis- 
tent to the more systematic study by Brune et al£ With 
a large corner diffusion constant, on the other hand, the 
scaling form of the island density changes drastically, be- 
cause small clusters of sizes less than 3 can migrate on 
the surface, i.e. i* = 3. The density is well fitted to the 
scaling n ~ F 3 ^ 7 . The increase of the critical island size 
i* is also reflected in the island size distribution: 2 -' 31,3 - 2 
The result is summarized in a scaling plot of the island 
density, and discussed in the last § 5. Even though the 
scaling of the island density to the deposition flux agrees 
to the rate equation analysis, the dependence on the dif- 
fusion constants fails to follow the expectation 



II. MONTE CARLO SIMULATION 

Our main aim is to explore the systematics and the 
universal features of the morphology change of islands 
on a singular surface as the diffusion mechanism changes, 
rather than the precise reproduction of some specific ex- 
perimental observations. For this purpose, we adopt the 
simplest model such that atoms are depositing on a (100) 
surface of a simple cubic lattice with a solid- on- solid 
(SOS) restriction. The SOS restriction is valid at low 
temperatures, common to MBE. 

Atoms are deposited randomly onto the substrate with 
a deposition rate of F monolayers per unit time. Here we 
choose the unit of time such that F = 1. The tempera- 
ture is assumed to be so low that all deposited atoms stick 
to the surface and never evaporate back into the ambient 
gas. An isolated adatom hops to one of its four near- 
est neighboring sites at a rate k s ; the waiting time of an 
adatom before the jump is (4fc s ) _1 , and the surface dif- 
fusion constant is given by D s = k s a 2 , where a is the lat- 
tice constant. As for the interlayer diffusion, which takes 
place in case a deposited atom lands on another adatom, 
no additional energy barrier is assumed for it to diffuse 
down. On the contrary, atoms at the edge are not allowed 
to hop up a layer. Therefore, Ehrlich-Schwoebel mound 
formation is supressed, and the layer-by-layer growth is 
expected. 

When an atom touches to other adatoms or islands, 
it makes a nearest neighbor bond to form a cluster. At 
very low temperatures, this cluster will never dissociate 
again. We assume that adatoms attach clusters irre- 
versibly, and no cluster dissociation takes place. This 
means that dimers are stable, and the critical island size 
is expected to be i* = 1. Nevertherless, if an edge atom is 
singly-bonded in the layer, i.e. if it makes only one bond 
with another adatom, it is possible to diffuse along the 
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(a) an ee atom 



(b) an ec atom (c) a cc atom 



FIG. 1: Diffusion processes of an adatom connected to an 
island with a single in-layer bond. k e is the rate of edge 
diffusion and k c is that of corner diffusion. 



island periphery at moderate temperatures. There are 
two possible diffusion jumps for a singly-bonded atom; 
an edge diffusion along the straight step edge, and a cor- 
ner diffusion or corner crossing round the outer corner. 
The rate of edge diffusion is k e and that of corner is 
k c , as is depicted in Fig^ Then, a singly-bonded atom 
can be classified by its possible motions in two direc- 
tions. It is classified as (a) an ee atom on a straight 
step if both movements are along the step edge, (b) an 
ec atom at the ridge of a corner if one is the edge dif- 
fusion and the other is the corner, and (c) a cc atom on 
a tip if both movements are across the corner, as shown 
in Fig^ The edge diffusion constant along the straight 
step edge is D e = k e a 2 , and we may define the corner dif- 
fusion constant as D c = k c a 2 . Precise values of diffusion 
constants depend on the energy barriers and the tem- 
perature. For edge and corner diffusion there may exist 
extra energy barriers in addition to that for the surface 
diffusion, and diffusion constants are probably ordered as 
D s > D e > D c > 0. If an edge atom is incoorporated 
into the kink site with more than two nearest neighbors, 
it is assumed to cease migration, for simplicity. 

Simulation starts from a clean substrate surface of a 
size L x L without adatoms. Hereafter we take the lat- 
tice constant a = 1. Therefore, the jump rates fc's and 
diffusion constants D's are the same. During the simula- 
tion, at some stage, there are No isolated adatoms, N ee 
edge atoms attached to a straight step edge, N ec atoms 
at the corner, and N cc atoms at the tip position (see 
Fig^l on the surface. These adatoms are sorted out in 
corresponding lists. The transition probability from this 
configuration to the next one is given as follows. The rate 
of the deposition is Pd = FL 2 , of the surface diffusion 
P s = 4A £> S , of the edge diffusion P e = (2N ee + N ec )D e 
and of the corner diffusion P c = (N ec + 2N CC )D C . The 
total probability of a state change in a unit of time is 
P t = Pd + P s + Pe + Pc- Thus, within a time interval 
dt = 1/Pt, one of the events takes place; the deposition 
with a probability Pddt, the surface diffusion with P s dt, 
the edge diffusion with P e dt and the corner diffusion with 
P c dt. For each diffusion process, the atom to be moved 
is picked up randomly from the list. After each state 
change, adatom lists have to be adjusted, for example, by 
eliminating those adatoms with more than two neighbor- 
ing bonds from the lists, or add an adatom that happens 
to have a single bond to the corresponding list. The time 
is increased by dt. Since the time unit is chosen such that 
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the deposition flux F = 1, the monolayer is covered at a 
time t = 1, in principle. Because of the stochastic fea- 
ture of the simulation algorithm, in practice, the actual 
time for the completion of a monolayer fluctuates around 
t = 1. 



III. ISLAND MORPHOLOGY 

Simulations are performed for systems of sizes 200 x 
200, and 1000 x 1000 under various combinations of dif- 
fusion constants, and some typical island morphology is 
shown in FigElat a coverage of 9 =0.1 monolayer (ML). 
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FIG. 2: Island morphology at various combinations of the 
surface, edge and corner diffusion constants, D s ,D e and D c . 
D s is fixed at D B /F — 10 9 . Parameters in parenthes repre- 
sents (D e /D a , D c /D s ). The system size is 200x200, and the 
coverage is 9 = 0.1 ML. 

From Figl^a) to (e) the edge and corner diffusion 
constants, D e and D c , are equally increased from to 
D s , by keeping the surface diffusion constant D s fixed 
at D s /F = 10 9 . Without peripheral diffusion, the is- 
land takes an irregular dendritic form or a fractal DLA- 
like shape, as shown in Fig^a). On increasing the 



edge and corner diffusion to D e = D c = 10~ 6 D S , den- 
drite branchs thicken and they extend mainly in the 
< 11 > direction, as shown in Fig[2Ib). Further in- 
crease of D e = D c = 10~ 5 D S let the islands compact 
in a hopper shape with a corner instability due to Berg 
effect^ as in Fig|^c). Since the corner is extending in 
the region with high monomer concentration, the comer 
captures many adatoms and the edge diffusion transfers 
them to the island center. But the edge diffusion is still 
weak and insufficient, and the growth of the central part 
is too slow to catch up the fast corner growth. Thus, 
the island grows mainly in the < 11 > direction, as in 
FigEIb) and (c). With sufficient edge and corner diffu- 
sion as D e = D c = 10~ 4 L\,, islands take square shape 
surrounded by rather straight {10} steps, as in Fig[2Id). 
So far island density seems independent on the island 
morphology. When all the diffusion constants are equal, 
as in Fig|2Ie), islands are square and very large, but their 
density decreases drastically. On this density change we 
shall discuss later in §4. The sequence of simulations 
from Fig[2Ib) to (d) almost corresponds to the param- 
eter range studied by Bales and Chrzanii by changing 
the temperature. The morphology variation is of course 
quite analogous to what they have obtained. 

We now study the effect of corner diffusion more pre- 
cisely, in a sequence of pictures shown from Figl^e) to 
(i). There, the edge diffusion constant is kept as large 
as D e = D s . and only the corner diffusion constant 
D c is decreased from D s to 0. In all these cases, is- 
lands are not irregular. With a large corner diffusion 
constant, D c > 10~ 4 D S , island are square with {10} 
steps, as shown in Fig(2[e) and (f). At D c = 10~ 5 D S 
in Fig[2Jg), they look round and rather isotropic, and for 
a still smaller D c = 10^ 6 D S in Figl^h) some islands are 
in diamond shape and some others are in a cross shape 
with thick arms. In the extreme case of Fig^i) without 
the corner diffusion, D c = 0, islands consist four sharp 
needles growing in the < 10 > direction. The same four- 
needle shape is observed previously by Zhong et ali& in 
a submonolayer growth and by Caspersen et al^ in a 
multilayer study. 

In a series of simulations Fig|2Ii)-(l) the corner dif- 
fusion is completely suppressed (D c = 0), and the ef- 
fect of the edge diffusion D e is exemplified. The series, 
in fact, continues back to Fig^a). With D e = D s in 
FigEIi) needles in cross-shape islands are rather stable 
and sidebranches are rare. On decreasing the edge dif- 
fusion to Dp = 10~ 4 L\, in Fig[21j), needles have many 
sidebranches growing in the < 10 > direction. Fur- 
ther decrement to D e = 1Q~ 5 D S leads to tip splitting in 
Fig.2(k). With still smaller D e = 10^ 6 D S clusters look 
like DLA aggregates. Comparison with Fig^b) and (1) 
reveals that the corner diffusion produces < 11 > den- 
drite in Fig[2fb), whereas without corner diffusion ag- 
gregates have a preference in the < 10 > direction in 
FigEIi)- Also without corner diffusion, the dendrite tip 
is pointed, whereas with a corner diffusion, the tip be- 
comes fat. At D e = D c = 0, the DLA looks isotropic. Of 
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course, it is known to be anisotropic asymptotically by 
reflecting the underlying square symmetry of the lattice, 
but the shape anisotropy appears only when the DLA 
grows very large. In the present simulation, islands coag- 
ulate to acomplish a layer growth before the asymptotic 
anisotropy appears. 

From these variety of island morphology the edge and 
corner diffusions are found to play different roles on the 
island symmetry. The edge diffusion smears out shot 
noise introduced by the irreversible attachment, and per- 
mits the formation of regular shape. The corner diffusion 
varies the orientational preference. Without the corner 
diffusion, those atoms on the tip positions or on the nar- 
row step edges cannot cross to the other sides of an island, 
and thus the supersaturation on the narrow {10} side in- 
creases. This leads to the enhancement in the growth 
velocity normal to the narrow {10} side, and due to the 
strongly anisotropic kinetic effect a needle tip is eventu- 
ally formed into the < 10 > direction. A corner diffusion 
allows the transfer of adatoms from one step edge to the 
other, and islands tend to be square. If the edge diffu- 
sion is weak, however, the adatom supply at the center of 
the step edge is meager, and the corner instability takes 
place. This leads to the dendrite growth in the diagonal 
< 11 > direction, as in FigEfb) and (c). 



IV. ISLAND DENSITY 
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FIG. 3: Double logarithmic plot of the island density n at 
the coverage 9 = 0.1ML versus the surface diffusion constant 
normalized by the deposition rate D a /F. The values of edge 
and corner diffusion constants, D e and D c , are (a) D e — D 3 
and D c = (no corner diffusion), (b) D e = D c = (with- 
out edge and corner diffusion), (c) D e = D c = 10" 4 D S , 
(d) D e = D c = W~ 3 D S , (e) D e = D c = l(n 2 Z) s , (f) 
D e = D c = 10 _1 D S , and (g) D e = D c = D s . In (a) the hori- 
zontal axis is displaced to the right by a factor 10, for a visual 
purpose. Lines represent scaling behaviors as n oc (D s / F)~ x . 
An exponent \ for straight lines is x = 1/3, and for dashed 
lines x = 3/7. 



In Fig [3 at a fixed surface diffusion constant D s , we 
observe that the island density is almost independent of 
the edge and corer diffusion constants, D e and D c . We 
said "almost", because in FigEfe) the island density is 
clearly different from the others. We now study the effect 
of peripheral diffusion on the island density. 

The rate equation analysis predicts that the island 
density n depends on the surface diffusion constant as 
in eqs.Q) and j2)l. 20 ' 21 i 24 ' 25 Experimentally, temperature 
dependence of the island density is used to estimate val- 
ues of various energy barriers by assuming the critical 
island size Simulations where the temperature is 

varied also has to assume i* to interpret their results. 22 
To study the relation of island morphology and its density 
by fixing i* = 1, we control directly on the rate of edge 
and corner diffusion, rather than tune many parameters 
as bond energies, energy barriers and the temperature. 

The island density n at the coverage 9 =0.1ML is plot- 
ted in Fig|2|as a function of the surface diffusion constant 
D s / F for various combinations of the edge and corner 
diffusion constants. The system sizes plotted are 200 2 
(open circles) and 1000 2 (filled circles). 

For the case without edge and corner diffusion, D e = 
D c = 0, islands are fractal as shown in Figl^a), at least 
at low densities with a large D s . The island density is 
shown in Fig|3fb). Data fits well with the scaling rela- 
tion n = 0.23(£> S /F)~ 1 / 3 . The result is consistent to 
that obtained by Venables et al& n = 0.25(D S /F ) _1/3 
from the rate equation analysis and that by Brune et al£ 



n = 0.27(L» S /F)- 1 027 / 3 for the DLA cluster in kMC sim- 
ulation at the coverage 0.12ML. For the DLA, in fact, 
another scaling law n oc [D s / F)^ 2 ^ i+D ^ is proposed, 
where Df = 1.71 is the fractal dimensioni2^i2Si2& From 
the present simulations of sizes up to 1000 2 , we cannot 
descriminate which of the two scaling laws is superior. 

For the case with a large edge diffusion D e = D s 
but without corner diffusion D c = 0, islands have cross 
shapes as in Fig^i). Even for this situation, the density 
n is compatible with a scaling law n = 0.25(D S / F)^ 1 ^ , 
as is shown in Fig^a). The plot is shifted by a factor 
10 to the right for the visual purpose. Otherwise, data 
overlaps almost completely with those for fractal islands, 
FigHb). 

With or without edge diffusion, the island density is 
found to follow the same scaling relation with the sur- 
face diffusion constant D s /F, obtained by a simple rate 
equation theoryi22i2i Inclusion of the corner diffusion, 
however, drastically changes the situation. By keeping 
D e = D c but at different ratio to the surface diffusion 
as (c) D e /D s = D c /D s = 10~ 4 , (d) 10" 3 , (c) 10" 2 , (f) 
10 _1 , and (g) 1, the relation between the island den- 
sity n and the surface diffusion constant D s /F varies 
as shown in Fig[3Ic)-(g). For (c) the island density at 
low D s /F < 10 10 is higher than that for (b) with frac- 
tal islands by about 20 percent, n = 0.28(D S /F )~ 1/3 . 
The density increase is due to the compactness of the is- 
lands, as shown in the surface configuration in FigQJb) 
for D e /D s = 10 -4 , compared to the fractal islands in 
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6 = 0.1ML. Ds/F = 10 9 and D e = D c . The ratio of the 
corner to surface diffusion constant D c /D s is (a) 0, (b) 10 -4 , 
(c) 1CT 3 , (d) 1(T 2 , (e) HT 1 , and (f) 1. 



FigEIa) for D c = 0. There are more space for adatoms 
to nucleate new compact islands. As for the D s depen- 
dence, the island density follows the same scaling law 
with an exponent 1/3, unless D s takes the largest value 
D s /F = 10 11 . There, the density becomes less than that 
of fractal islands, Fig|3fb). 

As D e and D c increase further from Fig|3fd) to (g), 
the island density starts to deviate from the curve 
FigHJc) at smaller D s /F. For the case Fig[3{d) with 
D e /D s = D c /D s = 10~ 3 , deviation takes place at about 
D s /F = 10 9 , for (c) at about D s /F = 10 7 . For (f) 
and (g) with larger D e = D c , the island density re- 
mains always less than that of fractal islands with the 
same D s . This tendency of density decrease is obvious 
by plotting the surface configurations in Fig[3| There, 
the coverage 9 = 0.1ML and the surface diffusion con- 
stant D s /F = 10 9 are the same, but the edge and corner 
diffusion constants, D e /D s = D c /D s , are varied. On in- 



creasing D c from Fig[3Jb) to (f), the number of islands 
decreases and the island size increases. 

From Fig|3fc) to (g) , the island density n seems to sat- 
isfy a new scaling relation to the surface diffusion con- 
stant D s /F as n oc (D s /F)~ x with an exponent \ ~ 0.43, 
when the density n is lower than that of fractal islands. 
In Fig|3]wc depict fitting lines with an exponent \ = 3/7, 
which is expected from the critical mobile island size 
i* = 3 and from cq.(0|. By fitting the island density 
in the scaling form with an exponent 3/7 as 



(3) 



the coefficient A depends on the peripheral diffusion con- 
stant D c = D e . By plotting A versus D c /D s in double 
logarithmic way as in Fig[S]the coefficient A is found to 
be proportional to (DJ D s )-° 195 . 




FIG. 5: The proportionality coefficient A as a function of 
the corner diffusion constant D c , which is equal to the edge 
diffusion constant D e . The line represents the curve A = 
0.38(D C /L> 3 )- 0195 . 
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FIG. 6: Island size distribution with (a) low and (b) high 
corner diffusion constant, D c . 

From Fig[U one further notices that as the corner dif- 
fusion increases, not only the total number of islands but 
also the island size distribution changes £L2L This ten- 
dency is apparent by plotting island size distribution in 
FigED Among the total Ni = nL 2 islands, there are 
N(s) islands of a size s. Since the average island size (s) 
changes as D c varies, the distribution of islands with a 
normalized size x = s/(s) has the relevant concern. Wc 
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define the probability P(x)Ax as the number of islands 
with sizes between x and x + Ax devided by the total 
number of islands Ni with Ax = 0.1. For each D c , size 
distribution of 10 samples are averaged. It is normalized 
as J °° dxP(x) = / °° dxxP(x) = 1, and corresponds to 
the scaled number density as N(s) = (6/ (s) 2 )P(s/ (s)) at 
a coverage At low L» c 's as D c = and 10~ 4 £> s in 

Fig|SJa), island sizes are widely spread with a broad peak 
such that there are very large islands with sizes about 2.5 
times the average size. This is the typical size distribu- 
tion for i* = On the contrary at high D c 's the 
island size distribution has a narrow peak around the av- 
erage size, and extends less than twice the average size in 
the region of large sizes. This is the typical size distribu- 
tion for i* > l i 27 i 31 i 32 When the clusters smaller then i* 
are dissociable, the island size distribution is proposed 3 ^ 
to have the form 



P l ,(x) = C l .x l exp(- 



(4) 



(i*ai») Qi *, and CV 



(i*Oj.) 



(i*+lK 



with 

r[(t* + 2) ai . 

r[(i* + i) 0i .] ~ v " ' ~ Oi.r[(«* + i)oi 

.(5) 

where Y means the Gamma function. The probabilities 
corresponding for i* = 1 and i* = 3 are drawn in FigjBfa) 
and (b) by continuous curves, respectively, and fit quite 
well with the simulation data, even though the dissocia- 
tion is not allowed there. In fact, the similar size distri- 
butions with mobile clusters are obtained in the previous 
Monte Carlo simulation^ There, however, no peripheral 
diffusion is granted, and the grown islands had fractal 
structures, in constrast to our square shaped clusters. 



V. DISCUSSIONS AND CONCLUSIONS 

We have seen that a large corner diffusion influences 
the island number density and the size distribution dras- 
tically. Without or weak corner diffusion, the island 
density varies in proportion to (D s /F)~ x with an expo- 
nent x =1/3, irrespective of the island morphology. This 
means that the nuclei are formed at the initial stage of 
deposition, as soon as the two adatoms collide with each 
other. They are immobile and act as a center of nu- 
cleation. Therefore, the randomness of the two-adatom 
encounter is frozen in. After sufficient density of nuclei 
is formed, further deposited adatoms are incorporated 
into the preexisting clusters. They develope into various 
shapes, depending on the edge and corner diffusion con- 
stants. Since the nuclcation takes place at random, the 
island sizes are distributed rather broad, reflecting the 
randomness in the island separation. In any case, this 
senario permits the island density to be independent of 
the morphology: The island density is determined long 
before the island shape appears. 

With a large corner diffusion, on the other hand, in- 
spections of the time evolution reveal that clusters of 



small sizes i < 3 = i* become mobile. Since atoms in 
a cluster smaller than the critical size i* = 3 can move 
around each other via the corner diffusion, small clusters 
migrate around randomly on the substrate surface. Dur- 
ing this cluster diffusion, they coagulate with isolated 
atoms and as a result the island density decreases. Of 
course, islands larger than the size i* =3 can migrate 
in some configurations, but when compact islands with 
doubly or more bonded atoms are formed, they cease 
migration on the substrate surface, and start to act as 
a nucleation center. During the further stage of cluster 
growth, the adsorbed atoms form small mobile clusters 
and are incorporated into islands in group. 

This situation is treated by Villain et ali24 in terms of 
the rate equation, and the island density n is expected 
to follow the scaling relation eq. J5J) . Since the trimcr can 
migrate on the surface, i* — 3, n oc F 3 / 7 is expected. 
This is the reason why we fitted the density with a scal- 
ing exponent \ = 3/7 in FigEl This explains the in- 
dependence of n, but as for the dependence on diffusion 
constants, we have to know D 2 and D 3 in terms of the 
present edge and corner diffusion constants, D e and D c . 

There is an exact calculation by Sanchez et al~ of 
diffusion constants of small clusters mediated by a pe- 
ripheral diffusion. The result is simple for a dimer and a 
trimer on a square lattice as 



Do = 



Da 



D P D r . 



3(D e + D c ) 



(6) 



For a dimer to diffuse, one of the adatoms which compose 
dimer has to cross round the corner. For a trimer, for 
example, in an extended form to move an edge atom has 
to cross the corner and to slide on the edge to another 
corner. Without edge diffusion, the center of mass of a 
trimer cannot move. In the case shown in FigO D e = D c 
and thus D% = D c /6. Inserting these exact results into 
the expression @ obtained by rate equation, the island 
number density is expected to follow the relation 



n oc 



F 3 
D S D 2 



1/7 



(7) 



If this scaling holds, the coefficient A in eq.@ should 
be proportional to A oc {D c / D,)- 2 ' 7 = (D c /D s )~ - 286 , 
instead the exponent -0.195 found in FigEl Therefore, 
the dependence of n on cluster diffusion constant does 
not seem to be properly given by the rate equation. 

To check the scaling behavior of the island density 
with the crossover from the monomer- to the trimer- 
limited form, let us study the scaled island density h = 
n(D s / F) 1 / 3 . It should be constant around 0.23 for small 
edge and corner diffusion constnat D e = D c . For large 
D e = D c , it should be proportional to (D 2 F/ D 3 ) 2 / 21 , 
according to the result of the rate equation analysis J3J. 
Therefore, we try the plot h versus F\ = D 2 F/D 3 , but 
the data does not collapse on a single universal curve. 

From our simulation result, the island density behaves 
as n = A(D S /F)~ 3 / 7 with an D c -dependent coefficient 
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FIG. 7: The scaled island density y = n(D s /F) 1/3 as a 
function of the scaled corner diffusion constant x = D 2 /D S F. 
The dashed curve represents small x asymptotic, y — 0.23 + 
0.05a? " 1 , and the continuous line the large x asymptotic form 
y = 0.38x~ 2/21 . 



A = 0.38(D c /D s )-° 195 for large D c , as depicted in 
FigEl Since the exponent 0.195 is approximately equal 
to 4/21, the density is expected to be scaled by the vari- 
able x = D^/FD S at a large D c more properly. The 
scaling plot is shown in Fig0for a wide range of combi- 
nations of diffusion constants as D s /F = 10 5 ~ 10 11 and 
D e = D c = 10~ la D s ~ D s , and one observes a nice data 
collapse on a universal curve. It is fitted as 



Ds 
F 



1/3 



0.23 + 0.05a; 01 
0.38a;- 2 / 21 



for x 
for x 





CO, 



(8) 



as shown in Fig[7| 

The effect of the edge and corner diffusion during the 
homoepitaxial nucleation growth on a singular surface is 
now summarized as follows. As for the island morphol- 
ogy, these peripheral diffusion makes the island regular 
and compact. With only the surface diffusion, the irre- 
versibly nucleated islands have a fractal shape charac- 
teristic to DLA. With the edge diffusion, islands become 
thick and compact. The corner diffusion affects the shape 



anisotropy. One may typically say that without the cor- 
ner diffusion dendritic needle pointing in the < 10 > di- 
rection appear, whereas with it squares with {10} faces 
appear. 

The density of islands is mainly determined before the 
shape appears. Without or with a weak corner diffusion, 
the critical island size is i* = 1, and the density n is re- 
lated to the surface diffusion constant D s normalized by 
the deposition rate F as n oc (Ds/F)' 1 / 3 , irrespective 
of island morphology. With a large corner diffusion con- 
stant, on the other hand, small clusters i < 3 can migrate 
randomly on the substrate surface, and another scaling 
relation n oc F 3 / 7 / D b J 21 D A C ' 21 = [F / D s f/ 7 (D s / D C ) A / 21 
is found appropriate. The flux independence agrees with 
the rate equation theory with i* = 3, but the dependence 
on diffusion constants D c contradicts the theory. The 
reason of this discrepancy is not clear. But, since the 
rate equation is a mean-field type approximation which 
neglects spatial correlation, it may be fortuitous that F 
dependence agrees in the simulation and the rate equa- 
tion theory. 

Recently, the mobility of large clusters is proposed to 
be very effective in selecting the island size and in ar- 
ranging islands in order in heterocpitaxial growth^ Our 
model allows the motion of only small clusters less than 
the size i* = 3 by the corner diffusion, and some size 
selection is observed. Due to the absence of any direct 
interaction between islands, however, the ordering in is- 
land positions cannot be achieved. Our model neither 
contains the size-limiting effect, and the coarsening takes 
place as the corner diffusion increases. 
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